Finite— momentum condensate of magnetic excitons in a bilayer quantum Hall system 
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We study the bilayer quantum Hall system at total filling factor Vt — 1 within a bosonization 
formalism which allows us to approximately treat the magnetic exciton as a boson. We show that 
in the region where the distance between the two layers is comparable to the magnetic length, 
the ground state of the system can be seen as a finite-momentum condensate of magnetic excitons 
provided that the excitation spectrum is gapped. We analyze the stability of such a phase within the 
Bogoliubov approximation firstly assuming that only one momentum Q is macroscopically occupied 
and later we consider the same situation for two modes ±Q. We find strong evidences that a first- 
order quantum phase transition at small interlayer separation takes place from a zero-momentum 
condensate phase, which corresponds to Halperin 111 state, to a finite-momentum condensate of 
magnetic excitons. 

PACS numbers: 73.21.Ac, 73.43.Cd, 73.43. Lp, 73.43.Nq 



I. INTRODUCTION 

A bilayer quantum Hall system (QHS) consists of two 
two-dimensional electron gases (layers) separated by a 
small distance d under an uniform magnetic field B per- 
pendicular to the layers. Among the several possible con- 
figurations, we consider the one where each layer has fill- 
ing factor v — n<f>o/B = 1/2, such that the total filling 
factor vt = 1/2 + 1/2 = 1. Here, n is the electronic 
density o f ea ch layer and <po — hc/e, the magnetic flux 
quantum P^i 

The system is characterized by two parameters: the 
ratios d/£ and Asas/^c- Here, £ — ^JhcjeB is the 
magnetic length, the characteristic length scale of QHSs, 
Asas is the electron interlayer tunneling energy, and 
E c = e 2 /e£ is the characteristic Coulomb energy with 
e being the dielectric constant of the host semiconduc- 
tor. Although Asas an d the distance d are fixed for a 
given sample, the ratio d/£ can be modified by changing 
the magnetic field B and then adjusting the electronic 
density in each layer in such a way that the configura- 
tion vt — 1/2 + 1/2 = 1 is restored. Interestingly, a series 
of measurementsSlEl has shown that for d < d c w 1.8 •£, 
the bilayer QHS behaves as a single-layer QHS at v = 1, 
while for d > d c , as two independent two-dimensional 
electron gases at v = 1/2. In spite of the fact that the 
experimental data indicate a continuous transition be- 
tween these two situations, the so-called incompressible- 
compressible quantum phase transition, from the theo- 
retical point of view it is not clear whether the system 
undergoes a second-order quantum phase transitiorP or 
a first-order one smeared out by disorder!^ 

The ground state of the bilayer QHS at vt = 1 is 
well understood in two limiting cases: for small d/£, it 
can be described by the (incompressible) Halperin 111 
wave function, 9 while in the very large d/£ region, by 
two in depen dent (compressible) composite fcrmion Fermi 
liquids GSUil Interestingly, the Halperin 111 state can be 



seen as a Bose-Einstcin condensate (BEC) of magnetic 
excitons, where the electron and the hole are in dif- 
ferent layers .121 This analogy motivated us to employ 
the bosonization scheme^ 3 -' to study the bilayer QHS at 
vt = 1. Our main finding in this first studjWwas that 
a zero-momentum BEC of magnetic excitons is stable 
only for d < 0A£ (zero interlayer tunneling case). Such 
a result is in quite good agreement with the exact diago- 
nalization calculations on finite size systems, which show 
that the overlap between the exact ground state and the 
111 state is close to unit only for d < 0.5 ^liSUSI 

Although much theoretical workplSliSBS] has been de- 
voted to the intermediate region, d ~ £, so far there is no 
consensus about the nature of the ground state. For in- 
stance, a (pseudospin) density wavepl mixed Fermi-Bose 
trial wave functions,^ and a (pseudospin) spiral state^ 
have been proposed as possible candidates. A proper 
description of the ground state in the intermediate d/£ 
region is important since it will help us to determine the 
nature of the incompressible-compressible phase transi- 
tion. 

In this paper, we re visit the bilayer QHS within the 
bosonization formalisn j 13 ! 14 ! focusing on the intermedi- 
ate d/£ region. We propose that within this bosonic 
scheme the ground state of the system can be seen as 
a /mite-momentum BEC of magnetic excitons. We show 
that this is indeed a possible phase of the effective bo- 
son model that we have derived in Ref. 1141 provided that 
the (neutral) quasiparticle excitation spectrum is gapped. 
Our results also indicate that the instability of the zero- 
momentum BEC of magnetic excitons at d = QA£ re- 
ported in Ref. [14] indeed corresponds to a first-order 
quantum phase transition from a zero-momentum BEC 
of magnetic excitons to a finite-momentum one. 

Our paper is organized as follows: In Sec. [Hj we in- 
troduce an interacting fcrmion model to describe the bi- 
layer QHS, summarize the bosonization method,^ and 
recall the main steps to derive the effective boson model 
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FIG. 1. (Color online) Schematic representations: (a) Bilayer QHS. B is the magnetic field and d the distance between the 
two layers, (b) Zero-momentum BEC of magnetic excitons and (c) finite-momentum BEC of magnetic excitons with |^Q| = 1. 
Only the lowest Landau levels corresponding to the t and 4- layers are shown, m is the guiding center quantum number which 
labels the degeneracy of each Landau level. 



from the original fermionic one. We also comment on 
the motivation for considering a finite-momentum BEC 
of magnetic excitons as the ground state of the bilayer. 
Sec. |IH| is devoted to the analysis of the effective bo- 
son model within the Bogoliubov approximation assum- 
ing that the ground state is given by a finite-momentum 
BEC of magnetic excitons where the momentum Q = Qx 
is macroscopically occupied. The ground state energy 
and the (neutral) quasiparticle excitation spectrum are 
calculated. Here, evidences that a first-order quantum 
phase transition takes place at small d/i are found. In 
Sec. |IV[ we perform a similar analysis but now consider- 
ing that two modes, ±Q with Q = Qx ^ 0, are macro- 
scopically occupied. We show that |ZQo|, the magnitude 
of the momentum associated with the lowest energy con- 
figuration, increases with d/l. Some additional features 
of a BEC of magnetic excitons are shown in Sec. |Vj In 



Sec. VI we compare our results with previous ones and 



comment on their consequences for the bilayer QHS at 
Vt = 1. A short summary with the main results closes 
the paper. The fact that density fluctuations can account 
for the definition of boson operators for the bilayer QHS, 
comparison with alternative bosonic schemes used to de- 
scribe the bilayer, and some details of the calculations 
can be found in the Appendices. 



II. MODEL 

Let us consider a two layer system composed of N 
electrons moving in the (x,y,z = 0) plane and N in 
the (x, y, z = d) plane under an external magnetic field 
B = Bz, Fig. flja), at zero temperature. We introduce a 
pseudospin index a =t, I in order to label each layer. We 
also assume that the B field is strong enough such that 
the electrons are fully spin polarized (frozen electronic 
spin degree of freedom) and that the Hilbert space of 
each layer is restricted to the corresponding lowest Lan- 
dau level. The configuration vt = v^+^i = 1/2+1/2 = 1 
is realized by setting the degeneracy of each Landau level 
iV = 2N. 

The Hamiltonian of the system has only two terms 
(since all electrons are restricted to the lowest Landau 
level, the kinetic energy is a constant and can be ne- 



glected): 



H = H 7 



Hr 



(1) 



Here, Ht describes the electron tunneling between the 
two layers, 



H 7 



AsAS^ C mt Cm ^ 



H.c. 



(2) 



and Hi is the Coulomb interaction term (we set the sys- 
tem area A = 1), 



>-\X 



Hl= 2^ ^ «a/j(*)pa(k)P/9(-k) (3) 

with k — |k|. Asas is the electron interlayer tunneling 
energy, c\ na creates an electron with guiding center m 
in the lowest Landau level of the a layer, Fig. Iljb), and 
/3a (k) is the Fourier transform of the a-electron density 
operator projected into the lowest Landau level, i.e.™ 



Pa(k) 



-(ik) 2 /4 



/ , k m]rn /(tk)c m 



(4) 



The function G m ^ m > (x) is defined in the Appendix C of 
Ref.HHl Finally, ' 



v n (k) = v u (k) = v A (k) = 
v n( k ) = v lt( k ) = v E(k) = 



2ne 2 

27TC 2 

IT 



,-("072 



-(lk) 2 /2 -kd 



(5) 



are, respectively, the Fourier transforms of the intralayer, 
va(t) = e 2 /er, and interlayer, ve{t) — e 2 /eVr 2 + d 2 , 
electron-electron interaction potentials with r = |r|. 

On can show that Eq. (J2J) can be written in terms of 
the 2;-component of the pseudospin density operator, i.e., 



Ht 



-A SAS S x (k = 0), 



(6) 



while Eq. pi , in terms of the total electron density op- 
erator p(k) = pf (k) + Aj.(k) and the z-component of the 
pseudospin density operator S z (k) = [pt(k) — /°i(k)]/2, 
namely 

H i = \ E Mk)p(k) P (-k) +2j2vz(k)S z (k)S z (-k), 

(7) 
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FIG. 2. (Color online) Schematic representations: (a) Quan- 
tum Hall ferromagnet |FM), the ground state of the single- 
layer QHS at v = 1, and (b) an electron-hole pair excitation 
(magnetic exciton) with momentum |/qj = 1 above |FM). 
Only the spin up and spin down lowest Landau levels are 
shown. E z is the Zeeman energy and m is the guiding cen- 
ter quantum number, (c) Semiclassical representation of an 
electron-hole pair in the symmetric gauge. Particles move 
along cyclotron orbits (solid blue circles) whose center are at 
one of the possible guiding centers (dashed red circles). The 
vectors R e and R^ correspond to the guiding center posi- 
tion of the electron (solid circle) and the hole (open circle), 
respectively. 



with 



v o/z(k) = -[v A {k)±v E {k)} 



-(lk) 2 /2 



ck 



\±e~ kd ) 



In the following, we focus on the zero tunneling case, i.e., 
we set Asas = which yields H — Hi. 



plicitly take into account the electronic spin. The ground 
state of the system, the so-called quantum Hall ferro- 
magnet |FM), is illustrated in Fig.[2la). It is possible to 
show that the neutral elementary excitations above this 
state, electron-hole pairs also known as magnetic exci- 
tons, Fig.pTb), can be approximately treated as bosons. 
More precisely, we can define the following bosonic oper- 
ators: 

where c\ na (c mCT ) is a creation (annihilation) operator 
for an electron in the lowest Landau level with guiding 
center m and spin a. The boson operators ([9| obey the 
canonical commutation relations [&L 6JJ = [6 q , 6k] = 
and [bq, fojj = 5 q k once some conditions are fulfilled. The 
state 6 q |FM) corresponds to a magnetic exciton with mo- 
mentum q, Fig. [2Kb). Within this framework, the elec- 
tron density operator and the z-component of the spin 
density operator read 

p(k) = 6 kt0 N* + 2*Y. sin ( k A q/ 2 ) & k+q 6 q' ( 10 ) 
q 

S,(k) = i<J fc ,oJV» - J2 cos ( k A q/ 2 ) fe k+q 6 q ( n ) 
q 

with k A q = £ 2 z ■ (k x q). 

It is easy to see that, in principle, the bosonization 
scheme outlined above can be employed to study the bi- 
layer QHS at Vt = 1> once the pseudospin a is identi- 
fied with the electronic spin quantum number a of the 
single-layer QHS at v = 1 [compare Figs UXb) and [2|a) 
and recall that we consider that the electrons are com- 
pletely spin polarized in the bilayer QHS]. Since the 
bosons b are defined with respect to a reference state, 
the quantum Hall ferromagnet |FM), the bilayer QHS at 
vt = 1/2 + 1/2 = 1 corresponds to a system with iV$/2 
bosons, as illustrated in Fig. [IJb). 



Bosonization formalism 



B. Effective boson model 



We study the interacting fermion model Q within the 
bosonization formalisrrpS that was recently developed by 
two of us among others. Although such a scheme was 
originally proposed for the single-layer QHS at v = 1, it 
is possible to show that it also holds for the bilayer QHS 
at vt = 1 , see Appendix |Aj We now briefly summarize 
the bosonization method and refer the reader to Ref. [T3] 
for more details. In Appendix |Bj we briefly comment on 
some alternative bosonic descriptions employed to study 
the bilayer QHS. 

Let us consider the single-layer QHS at v = 1. We re- 
strict the Hilbert space to the lowest Landau level and ex- 



Let us now follow the lines of Ref. [13] and map the 
original interacting fermion model ufh into an effective 
interacting boson model. Substituting Eqs. ( 10 ) and ( 11 ) 
into Eq. (|7| and normal ordering the result, we arrive at 



H B =Y,"qbqbq+ ^ V k ( P , q)& k+p & q _ k & q 6 p . (12) 
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FIG. 3. (Color online) Dispersion relation of the free bosons 
(in units of e 2 /0, Eq. |l3j, for d/t 
(from bottom to top at £q = 0). 



0, 0.5, 1, 1.5, and 2 



is the dispersion relation of the free bosons (see Fig. [3]), 
with Jq(x) denoting the Bessel function of the first kind 
and 

^k(p, q) = 2v (k) sin(k A p/2) sin(k A q/2) 

+ 2i^(fc)cos(kAp/2)cos(kAq/2) (14) 

is the boson-boson interaction potential. In the follow- 
ing, instead of Hg , we consider 



K = H B - f iN, 



(15) 



which explicitly includes the chemical potential jjl. Here, 
N = J2q ^q^q ^ s tne number operator for bosons. 



C. Finite— momentum BEC of magnetic excitons 



In Ref. HH we analyzed the interacting boson model 



( 12 1 assuming that the bosons b condense in their lowest 
energy state, the q = mode, and showed that such a 
state is stable only for d < 0A£. The good agreement 
between our results and exact diagonalization calcula- 
tions, see the Introduction section, tells us that not only 
is this zero-momentum BEC a good approximation for 
Halperin 111 state, but also that the bosonic formalism^ 
is indeed quite appropriate to describe the bilayer QHS 
at v<r = 1. Therefore, it might be possible to describe 
the decrease of the correlations between the two layers 
as d/£ increases, i.e., the intermediate d/t region, using 
solely the bosonic degrees of freedom. In this case, what 
should be the form of the ground state in terms of the 
bosons b for d ~ 11 

In order to construct the new ground state, we should 
recall some properties of the magnetic excitons. As men- 
tioned above, the state 6_|FM) corresponds to a magnetic 
exciton with momentum q, which is nothing but a suit- 
able linear combination of electron-hole pairs above the 
|FM) state, see Eq. (J9l). The momentum q is canonically 
conjugate to the vector R = (R e + R/j)/2p3 where the 
vectors R e and R^ denote respectively the position of 



the guiding centers of the electron and the hole as illus- 
trated in Fig. [2] (c). Interestingly, it is also possible to 
show that [see Eq. (2.16) in Ref.lMj 



(FM|& q (R e - R h )6 q |FM) 



qx z, 



(16) 



i.e., the (relative) distance between the guiding centers of 
the electron and the hole which constitute the magnetic 
exciton is ex q. Note that this is an unusual relation 
between momentum and distance. Therefore, a boson b 
with q = can be seen as an electron-hole pair both 
localized in the same guiding center, while for a boson 
b with q^ 0, the electron and the hole are in different 
guiding centers. 

A zero-momentum BEC of magnetic excitons is then 
characterized by a large number of (interlayer) electron- 
hole pairs where each electron is very close (in the guiding 
center sense) to its partner hole as depicted in Fig. [I] (b) . 
Since this is the smallest distance between the electron 
and the hole, such a feature indicates that the two layers 
are highly correlated, corroborating the relation between 
the zero-momentum BEC and the 111 state. Therefore, 
in order to decrease the coupling between the two layers, 
we should, in principle, consider a state constituted of a 
large number of electron-hole pairs where now each elec- 
tron is a little bit displaced from its partner hole. This 
situation is nothing but a finite-momentum BEC, where 
the bosons macroscopically occupy a finite Q mode, for 
instance, the one with tQ = |^Q| = 1, Fig. [11(c). Given 
such a relation between the momentum Q and interlayer 
coupling, we also expect that the larger tQ, the lower the 
correlation between the two layers. 

These are the key points which motivated us to pro- 
pose a finite-momentum BEC of magnetic excitons as a 
possible ground state for the bilayer QHS in the interme- 
diate d/t region. In the next two sections, we study the 
stability of this state at two different levels of approxi- 
mation. 

As a final remark, we should note that although finite- 
momentum BECs have been recently discussed in the 
context of ultracold Bose gases (see, for instance, Refs.l25l 
and I26p . our motivation to consider such a phase is 
mainly due to the properties of the magnetic exciton as 
explained above. 



III. ONE-MODE APPROXIMATION 



In this section, we analyze the effective interact- 



ing boson model (12) within the so-called Bogoliubov 
approximation 27 assuming that the ground state is given 
by a finite-momentum BEC with the Q = Qx mode 
macroscopically occupied. We hereafter refer to this pro- 
cedure as the one-mode approximation. Although the Q 
mode is not the lowest energy single-particle boson state, 
see Fig. [3j we show that such a BEC is indeed a stable 
solution for certain values of d/t provided that the excita- 
tion spectrum is gapped. Here the boson-boson interac- 
tion potential ( 14 ) plays an important role in the stability 



of this phase. In the following, we consider 0.1 < IQ < 2 
and 0.1 £< d<4£. 

Before continuing, some words about the approxima- 
tion scheme are here in order: since the single-particle 



boson energy (13) has cylindrical symmetry, w q = uj qi 
there is no reason to expect that the bosons will condense 
in only one particular momentum Q = Qx. In principle, 
the bosons b could even condense in more than one mode 
q as long as q = Q. However, such an approximation 
is the simplest one which allows us to verify whether a 
finite-momentum BEC of magnetic excitons is indeed a 
stable phase via quite accurate and well controlled cal- 
culations. This is the idea of the procedure adopted in 
this section and in the next one. Later, in Sec. |VI[ we 
will discuss which of the features found here could be 
displayed by the bilayer QHS and also how the results 
derived from these two initial considerations could guide 
us to propose a more elaborated approximation scheme 
to study such a finite-momentum BEC. 

We start by replacing 6L = 6q = (bL) = (6q) — >• V^o 
in Eq. (15), where No is the (macroscopic) number of 
bosons in the Q mode. Keeping only terms with two 
bosonic operators, one finds after some algebra that 



le+&J 



K=K +-J2 [tq«Q + q «Q -q 

q#0 



e~b 



»,t 



q "Q-q"Q-q 



v qV"Q+q"Q-q 



^Q-q^Q+q^ 



(17) 



where 



K Q = (cjq - n)N a - - Y^ e 



2 ^-^ ~ q : 

q#0 



£ q — W Q±q — M+ ^q; 



A q = 2JV U q (Q,Q). 



(18) 



The quadratic Hamiltonian (17) can be diagonalized with 



the aid of the canonical Bogoliubov transformation 



bn^„ = U„aU^„ - VqOQ-q, 



3 Q+q 
fe Q- q 



x q"Q+q 

U q a Q-q 



UqOQ+q, 



(19) 



which yields 
K = K 



Here 



»+jE (°q a Q+q a Q+q + °q a Q-q a Q-q 

q#0 

Ko + ^ ^a^aq. 
q^Q 



n± = ±A q + n q , 



(20) 



^q — y e q - A q , 

1 / _u 

2 



A q = 9 (4 - £ q) 



K = (wq - n)N --^(e q - Q q ) , 



(21) 



the quasiparticle dispersion relation is given by 

"q — "q-Q' 

and the Bogoliubov coefficients obey 



(22) 



1 

2 + 

A 

2fi c 



2fi q ' 



c q 

2n n 



(23) 



The chemical potential fi can be obtained from the sad- 
dle point condition dK /dN = 0: since dcq/dN = 
dXq/dNo = Xq/No, one can show that 



/' 



W Q + 7^H VqK 



N 



WQ + Ho- 



(24) 



q#0 



No follows from the conservation (on average) of the total 
number of bosons N B = £ q (& q & q ) = N^/2 = l/4ir£ 2 : 



bosons in the condensate is 



from Eqs. (|19|), one finds that the relative number of 

(25) 



n 



N () 



N f 



i-E" 



q#0 



Finally, the ground state energy E (Q,d) = K + fi{N) 
reads 



E Q (Q,d) _ K 



N B 



N b 



Jq+Mo = WQ+^o(l-no)-/oi (26) 



With Ioi = gA- E„^0 ( 



^q^O V c q 



fir 



Once jj, and no are 



known for fixed Q and d/£, the quasiparticle spectrum 
fi q and the ground state energy ( 26 ) are completely de- 
termined. 



A. Zero-momentum BEC 

Before proceeding, we would like to briefly recall the 
results from our first analysis of the effective boson model 



( 12 ) reported in Ref. [TJ] 

By setting Q = and ^ = in the above equations, 
we recover Eqs. (8) and (9) of Ref. [HI The choice no = 0, 
based on the one-loop approximation, 28 yields a gapless 
excitation spectrum for the zero-momentum BEC phase 
[see Fig. 2(a) from Ref. 1 14). in agreement with the Gold- 
stone theorem. 



q#0 



We also find that the ground state energy (26) in- 
creases with d/£, Fig. H (long dashed black lineJT and 
that the relative number of bosons in the condensate 
no, Eq. (25), decreases rather fast as d/l increases, 
Fig. [5](a). Indeed, such a behavior of n led us to in- 
clude into the description the quartic terms in boson op- 
erators of the Hamiltonian ( 12 ) neglected in the Bogoli- 
ubov approximation.— Considering these quartic terms 
in the so-called Popov approximation,— we showed that 
the self-consistent equations admit solutions only for 
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FIG. 4. (Color online) Ground state energy per boson (in units of e 2 /t£) as a function of d/£: (a) one-mode approximation and 
(b) two-mode approximation. £Q — 0.5 (dotted magenta line), 1 (solid red line), 1.5 (dashed blue line), and 2.0 (dot-dashed 
green line). Long dashed black line: Q — 0, one-mode approximation with po = 0, see text for details. Inset: details of the 
corresponding main plots focusing on the large d/l region. 



d < d c o = 0A£. Here, we revisited the problem and 
perform more accurate numerical calculations. We find 
that d c o = 0.56 I, w hich i s even closer to the exact diago- 
nalization estimate d 15 * 16 ! mentioned in the Introduction. 



B. Finite— momentum BEC 

Let us now consider £Q ^ and discuss numerical 
solutions of Eqs. (24)— (25 ). It is possible to solve the 

in the con- 



self-consistent problem for all values of 
sidered range as long as a finite (self-consistently cal- 
culated) fi is allowed and d is larger than a minimum 
value d m in. This feature is exemplified in Fig.Wla), where 
we show the ground state energy ( 26 1 as a function of 
d/£ for IQ — 0.5, 1.0, 1.5, and 2.0. One can see that 
^min = 0.1, 0.3, 0.8, and 1.9 £ respectively for d = 0.5, 
1.0, 1.5, and 2.0 £. Note that the four configurations lie 
quite close in energy as d/t increases, but the ground 
state energy curves never cross each other. This behav- 
ior is also observed for all intermediate £Q values (not 
shown here), i.e., Eo(Q,d) increases with £Q for a fixed 
d/l. It is clear that a finite-momentum BEC is lower 
in energy than the zero-momentum BEC discussed in 
the previous section for d > 1.0 £ Interestingly, the 
Eq(Q — 0, d) and E (Q ^ 0, d) curves cross at a (small) 
critical layer separation d c \, indicating that a first-order 
quantum phase transition from a zero-momentum BEC 
to a finite-momentum one takes place at this critical 
value. Note that for configurations with 0.5 < £Q < 1.0, 
d c i is within the range 0.45 £ - 0.7 £, which includes the 
(updated) d c0 previously determined within the Popov 
approximation in Ref. 1141 

One important consequence of a finite /xo is that the 
dispersion relation of the (neutral) quasiparticles is now 
gapped. For instance, in Fig. |6l we show the excitation 
spectrum (22 1 along some particular momentum direc- 



tions for the configuration with £Q — 1 a,t d — 1.2£. The 
minimum gap A is at a momentum qA = —Qax, i-e., the 
angle between qA and Q is equal to n. For a fixed IQ, 
qa continuously increases with d/£. We also find that, 
for a given IQ, the gap increases with d/l as shown in 
Fig. [7] (dashed lines) . The fact that a gap opens up at d c \ 
provides further support for a first-order quantum phase 
transition at this critical layer separation. Finally, note 
that J7 q has no longer cylindrical symmetry, n q ^ £l q , 
which differs from the excitation spectrum of the zero- 
momentum BEC [Fig. 2(a), Ref. IT4] . This aspect and the 
peak in f2 q at q = Q are artifacts of the oversimplified 
one-mode approximation. 

In order to understand the behavior of the excitation 
spectrum at small momentum q, we should look at the 
nature of the elementary excitations. Recall that a bo- 
son b has an internal structure since it corresponds to 
an electron-hole pair. An elementary excitation of the 
magnetic exciton BEC can be seen as an electron-hole 
pair with momentum Q which is taken from the con- 
densate, broken and recombined again in a electron-hole 
pair but now with a momentum q 7^ Q. Apart from the 
corrections due to the boson-boson interaction potential, 



Eq. ( 14 1, the excitation energy f2 q is related to the differ- 



ence AEb between the binding energies of the pairs with 
momentum Q and q, namely 



o. r 



AE h = 



; Q 



where w q is the dispersion relation of the free bosons, 
Eq. (13). Let us firstly consider the zero-momentum 
BEC. In this case, the bosons are condensed in the low- 
est single-particle energy mode, Q = 0, and therefore 
lim q _j. (wQ — o>q) = which yields a gapless excitation 
spectrum, i.e., the system displays a Goldstone mode. On 
the other hand, in a finite-momentum BEC, the bosons 
are not condensed in the lowest single-particle energy 
mode. This is an important feature which implies that 
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FIG. 5. (Color online) Relative number of bosons in the con- 
densate no , Eqs. (251 and ( 42 1 , as a function of d/l: (a) zero- 
momentum BEC with fig = 0, see text for details, and finite- 
momentum BEC with IQ = 1.0 (solid red line), 1.5 (dashed 
blue line), and 2.0 (dot-dashed green line) within the (b) one- 
mode approximation and (c) two-mode approximation. 



lirriq^o^Q — W q ) 7^ 0, i.e., the Goldstone mode disap- 
pears. In other words, the internal structure of the boson 
b combined with a macroscopic occupation of a higher en- 
ergy single-particle mode leads to the disappearance of 
the Goldstone mode. Such a behavior reminds us of the 
excitation spectrum of a BCS superconductorPS We will 
return to this issue in Sec. IVI Al 

Finally, we find that the relative number of bosons 
in the condensate no, Eq. (25), is roughly independent 
of d/l and close to one. Such an aspect, illustrated in 



Fig. Ojjb), is related to the existence of a finite excitation 
gap which reduces quantum fluctuation effects compared 
with a gapless case (the zero-momentum BEC). The fact 
that no ~ 1 tell us that the Bogoliubov approximation 
is indeed quite reasonable to study a finite-momentum 
BEC phase, in contrast with the zero-momentum BEC, 
which requires a more involved approximation. 



IV. TWO-MODE APPROXIMATION 

So far we have considered that the bosons b condense 
in just one particular single-particle mode q = Qx. 
As mentioned in the previous section, since the single- 
particle boson dispersion relation ( 13 ) has cylindrical 
the bosons b could, in principle, 



symmetry, u; q = us q 
condense in more than one mode q provided that q = Q. 
In this section, we discuss such a possibility, in partic- 
ular, we assume that the BEC is split into two pieces: 
both q = ±Q modes, with Q = Qx and IQ ^ 0, are 
now macroscopically occupied. Again, the Bogoliubov 
approximation is employed to analyze the effective boson 
model (12). We hereafter denote such scheme two-mode 



approximation. 




FIG. 6. Dispersion relation of the (neutral) quasiparticles, 
Eq. p2J), (in units of e 2 /e£) for a finite-momentum BEC with 
£Q — 1 at d = 1.2^ along some particular momentum direc- 
tions within the one-mode approximation. 



by performing the substitutions 



(27) 



in Eq. (15). The equivalent of Eq. (17) is now given 
by Eq. (CI), see Appendix [Cj In ord( 



: er to diagonalize 
the Hamiltonian (CI), it is useful to introduce the four 
component vector 



q 



(« 



/,t 



Q+ q b -Q+ q & -Q-q & Q-q 



) 



Eq. (CI I can then be expressed in matrix form: 



q- 



where the 4x4 matrix 7? q reads 



H„ 



S q i q 



(28) 



(29) 



(30) 



with the 2x2 matrices A q and B q given by 



A q = ' Cq 



A, 



7q 



7q 



A r 



..!.(! B a = I ^ A<1 



Aq £q 



Here, we assume that both condensates have the same 
number of bosons and set Nq — No. The coefficients e q 



and A q are defined in Eq. (21) whil e 7 q , £ q , and A q are 
shown in the Appendix [C] sec Eqs. (C2) and (C3). 
The diagonalization of the 4x4 problem ( 29 ) is more 



involved than the 2x2 one corresponding to Eq. (17). 



Therefore, it is more convenient here to use the procedure 
described in Ref. [31] Since we are dealing with a bosonic 
system, instead of -ff q , one should diagonalize 



Here, we basically follow the lines of Sec. Ill and start 



IbH^ 



with 



In = 







(31) 




FIG. 7. (Color online) Minimum gap energy A (in units of 
e 2 /e£) of the (neutral) quasiparticle excitations as a function 
of d/t for a finite-momentum BEC with IQ = 1, 1.5, and 
2. Dashed lines: one-mode approximation; solid lines: two- 
mode approximation. 



The (positive) eigenvalues of the matrix (31 1 are 
where 



(32) 



C' 



D„ 



4 + A£ + 7 £-A£-£, 

A q (4 - A q ) + ( 7q£q - A q £ q ) 2 
Eq. ( 29 ) then acquires the form 



1/2 



where the 4x4 matrix H' reads 



//■■ [ ^q 
q I L 



with 



q- 



0+ 
fi: 



(33) 



(34) 



'q / \ " q 

and the new four component vector $ q is given by 



(35) 



(36) 



The relation between the two set of bosonic operators 
a±Q± q and 6 ±Q±q is 

* q = M q * q , where M q = ( fy M (37) 
with U„ and V„ being 2x2 matrices, 



U a 



_ fwi(q) u 3 (q) 
«2(q) w 4 (q) 



V a 



_ f Vi(q) w 3 (q) 
«2(q) u 4 (q) 



whose elements are the Bogoliubov coefficients. The com- 
plete expressions of the Bogoliubov coefficients Ui (q) and 



«i(q) are quite long and they can be found in the Ap- 
pendix JO 

Eq. ( 34 ) can be rewritten as 



K = K + ^2 ^q « q a q> 

q#±Q 
where the quasiparticle energy il q reads 

1 



<X 



fi 



Q+q 







-Q+q + ^Q+q 



n 



-Q+q 



(38) 



(39) 



and 



K = 2A 2 U2Q (Q, Q) + 2N (u Q - M ) 

+ iE("q+^q-<- e q)- 



(40) 



Again, from the saddle point condition OKq/ONo = 0, 
the chemical potential \x can be calculated: since from 
Eqs (fC2l) with A = A we have 



9e±/9^ = dXJdN = A q /A , dAJdNo = 0, 

d^/dN = 7q /A , d£JdN = C q /Ao , 
after some algebra, we find that 

M = Wq + Mo + 2N v 2 q(Q, Q) = wq + Mo + Mi- ( 41 ) 



The quantity moj see Eq. (C4|, is different from the one 



mode approximation expression, Eq. ( |24[ ) . From the con- 
servation (on average) of the total number of bosons 
n b = E q (^q) = 2A + £ q# ± Q (& q 6 q >, it follows that 
the relative number of bosons in the condensate uq is 
given by 



n 



Aq_ 



lEE^(q)- 



(42) 



i — 1 q 



Finally, it is easy to see that the ground state energy 
E Q (Q, d) = K + fj,(N) reads 



E (Q,d) K 



N F 



N f 



JQ + Mo + Mi 



^q +Mo(l - 2n ) + Mi(l - n o) - Io2 (43) 



with I< 



112 



Ec 



n- 



ft ) . Similarly to 



4N B ^q ( £ q + C q ""q , 

the one-mode approximation, Sec. |HI| we numerically 
solve the self-consistent Eqs. (41 ) and (42 ) and determine 



no and m fo r fixed ^Q and d/t. 

In Fig. J4^b) , we plot the ground state energy ( 43 1 as a 
function of d/t for three different configurations, the ones 
with tQ — 1.0, 1.5, and 2.0. Likewise the one-mode ap- 
proximation, the self-consistent equations can be solved 
only for d larger than a minimum value d m - ln . However, 
we now have <i m ; n = 0.6, 0.6, and 0.7 1 respectively for 
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FIG. 8. Magnitude of the momentum £Qo corresponding to 
the lowest energy configuration for a fixed d/l, two-mode 
approximation. Dashed line is a guide to the eyes. 



IQ = 1.0, 1.5, and 2.0, which differ from the one-mode 
approximation results. Comparing the ground state en- 
ergies obtained with both one-mode and two-mode ap- 
proximations for a given £Q, we clearly see that the latter 
is lower than the former: since w q = u q , macroscopic oc- 
cupation of both ±Q modes are equally likely. The sys- 
tem then profits from this fact by splitting the condensate 
into n = 2 equal pieces, binding them and lowering the 
total energy. Again, a finite-momentum BEC is more fa- 
vorable than a zero-momentum one for d > 1.0 £ and the 
Eq{Q i^ 0,d) and the E (Q = 0,d) [one-mode approxi- 
mation with /io = 0] curves cross at small d. In particu- 
lar, E (IQ — 1, d) and E (Q = 0,d) cross at d a \ « 0.681, 
in good agreement with the result derived in the previous 
section. Therefore, both one-mode and two-mode ap- 
proximations indicate that a first-order quantum phase 
transition may occur at small d/l. 

Concerning the large d region, we again find that the 
configurations with IQ / are quite close in energy 
but now, differently from the one-mode approximation, 
the different Eq(Q ^ 0,d) curves cross each other. For 
instance, Eq(£Q — l,d) and Eq(£Q — 1.5, d) cross at 
d « 3.1 £, see inset Fig. [41(b). Indeed, we find sev- 
eral crossings between the different ground state energy 
curves for 1.0 < IQ < 2.0 and 0.7 £ < d < 4£. In particu- 
lar, the magnitude of the momentum |£Qo| corresponding 
to the lowest energy configuration for a given d/£ is shown 
in Fig. [8] The fact that £Qo increases with d/£ corrobo- 
rates the scenario proposed in Sec. II C that the larger £Q, 
the lower the correlation between the two layers. Note 
that the one-mode approximation is not enough to cap- 
ture such a behavior. Moreover, the results also indicate 
that another first-order quantum phase transition may 
take place at larger d, from one finite-momentum BEC 
with small £Q Q to another one with a larger £Qq. In par- 
ticular, the transition £Q = 1.0 — ¥ £Qo = 1.1 occurs 
at d C 2 — 1.6 £, which is very close to the critical layer 
separation where the incompressible-compressible phase 
transition is experimentally observed. 3 Finally, we should 
mention that solutions for the £Q = 0.8 and 0.9 config- 
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FIG. 9. (Color online) Contour plot of the quasiparticle dis- 
persion relation f2 q , Eq. (391, (in units of e 2 /e£) for a finite- 



momentum BEC with IQ — 1 at d — 1.5 £, two-mode approx- 
imation. 



urations are also possible but, since they are very close 
in energy to the £Q = 1.0 configuration, we decided to 
neglected them in the above discussion. 

Differently from the one-mode approximation, here 
the inversion symmetry of the excitation spectrum f2 q , 
Eq. ( 39 ) , is preserved as exemplified in Fig. [9] for the 
finite-momentum BEC with £Q = 1.0 at d = 1.5£. 
Again, a finite (self-consistently determined) //o leads 
to a gapped excitation spectrum. Note that the mini- 
mum gap, which is larger than the corresponding one de- 
termined within the one-mode approximation, increases 
with d/£ for a fixed £Q, Fig. [7] (solid lines). Interest- 
ingly, for small d/£, the minimum gap is located at the 
origin (qA = 0) but, as d/£ increases, the position of the 
minimum gap abruptly changes to qA = IaV (the an- 
gle between qA and Q is now 7r/2, in contrast with the 
one-mode approximation result) and then £q& continu- 
ously increases with d/£. Such a behavior is exemplified 
in Figs. 10 for the £Q = 1.0 (upper row) and 1.5 (lower 
row) configurations. Also, the kinks observed in Fig. [7] 
(solid lines) are signatures of this abrupt change in qA ■ 
Finally, some words about the singularities of the exci- 
tation spectrum are here in order: we believe that the 
peaks in f2 q at ±Q (also found in the one-mode approx- 
imation) might be an artifact of the two-mode approxi- 
mation and that they may disappear as we increases the 
number of components n (even) of the finite-momentum 
condensate. We will return to this point in Sec. |VIB| 

Concerning the relative number of bosons in the two 
condensate pieces 2n , Fig. p)Fc), we can see that its be- 
havior is similar to the one found in the previous section: 
it is almost independent of d/£ and 2no is close to one. 
The latter indicates that the Bogoliubov approximation 
is indeed appropriate to study the finite-momentum BEC 
even if the condensate is split into more than one piece. 

As a final remark, we would like to comment on the 
fact that the results found in this section seem to con- 
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FIG. 10. (Color online) Quasiparticle dispersion relation fi q , 

directions for several values of the ratio d/£. Finite-momentum BEC with 



Eq. (39 1, (in units of e 2 /e£) along some particular momentum 
= 1.0 (upper row) and £Q = 1.5 (lower row), 



two-mode approximation [d/l increases from bottom to top at £q = 0). 



tradict Nozieresp^J who argued that a fragmentation of 
the condensate into two pieces costs a macroscopic ex- 
tensive exchange energy and therefore it is not favorable. 
A careful analysis shows that there is no contradiction. 
Let us denote E i and E 02 the ground state energy re- 
spectively obtained within the one-mode [Eq. (26|] and 
two-mode [Eq. (43 1] approximations. Using the fact that 
no [one-mode approximation, Eq. (25 ) ] is roughly equal 
to 2no [two-mode approximation, Eq. (42)], we have 



En 



En 



/xi(l - n ) - I02 + It 



01- 



Comparing the above equation with Eq. (4) from Ref. [32j 
we realize that Nozieres considerations only take into ac- 
count the first term in the above equation and completely 
neglect the other terms which, as we have seen, provide 
important corrections. In particular, for the bilayer QHS, 
symmetry considerations also indicate that a BEC split 
into two parts has lower energy than a single condensate: 
recall that the excitation spectrum obtained within the 
two-mode approximation is more symmetrical than the 
one derived in the one-mode approximation. 



PROPERTIES OF A BEC OF MAGNETIC 
EXCITONS 



It is easy to see from Eq. ( 10 1 that 



<p(k)> = <J fcj0 JV«, (44) 

which is valid for both one and two-mode approximations 



ciple, the BEC of magnetic excitons is an homogeneous 
phase (see discussion at the end of the section). Con- 
cerning the expectation value of the z-component of the 



pseudospin density operator (111, one shows (two-mode 
approximation) 

(S,(k)) = v^VVo exp(-KQ| 2 /4) (<5 k ,-2Q + S k , 2Q ) . 

.(45) 
Note that ( 45 ) vanishes within the one-mode approxima- 
tion since iVo — >• 0. 

Further insight into a BEC of magnetic excitons can 
be obtained by looking at the pair correlation function, 
which is defined aa^ 



w-i^xy^'Pfa)- 1 ]' 



where the static structure factor is given by 



5(q) 



1 
N 



(p(-q)p(q)) - NS qfi 



(46) 



(47) 



regardless the value of IQ. Eq. (44) implies that, in prin 



with p(q) being the Fourier transform of the electron 
density operator. The pair correlation function basically 
tell us the probability of finding an electron at the po- 
sition r giving that there is another one at the origin. 
The analytical expression of g(r), both at the one and 
two-mode approximations, is quite lengthy and can be 
found in the Appendix |T)| Here, we just comment on its 
numerical evaluation. 

In Fig. |11[ we plot the pair correlation function along 
one particular r-direction for the zero-momentum BEC 
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FIG. 11. Pair correlation functio n (|46[ ) along one particular r-direction. Dashed black line: zero-momentum BEC, one-mode 
approximation with /io = 0, Eq. ( |D5[ ). Finite-momentum BEC with (a) £Q = 1.0 and (b) £Q = 2.0, two-mode approximation. 
d — 1 (solide red line), 1.5 (long dashed blue line), and 2£ (dot-dashed green line). Insets: details of the corresponding main 
plots showing the behavior of the pair correlation function in the large r/l region. 



(one-mode approximation with /j,q = 0) and for the 
finite-momentum BECs with £Q = 1 and 2 at d = 1.0, 
1.5, and 2.0 £, calculated in the two-mode approximation. 
Within the approximations considered here, g(r) is d/£ 
independent for the zero-momentum BEC, see dashed 
line in Fig. [TT] In this case, the pair correlation function 
vanishes as r — > 0, indicating the existence of a correla- 
tion hole around the electron, and it is constant at large 
r, the same features displayed by the single-layer QHS 
at v = lp31 A distinct behavior is found for the finite- 
momentum BEC. Note that now g(0) ^ 0, indicating 
that two electrons (with different pseudospins quantum 
numbers) can be very close to each other, corroborating 
somehow the schematic picture for a finite-momentum 
BEC depicted in Fig. flic). Moreover, small oscillations 
at large r/£ are observed, which are characteristic of a 
composite fermion Fermi-liquid.^ These findings sup- 
port the proposal that the ground state of the bilayer 
QHS in the intermediate d/£ region can be described by 
a finite-momentum BEC of bosons b. 

As a final remark we would like to mention that in 
Ref. [35] it is shown that an exciton condensate has diag- 
onal long-range order. Interestingly, the average value 
of the density operator is constant and only the density- 
density correlation function has Fourier components of 
the type exp[— iK • (ri — r2)]. Therefore, based only 
on Eq. ( 44 ) , we would expect that the finite-momentum 



BEC of magnetic excitons corresponds to an inhomoge- 
neous phase. However, since the pair correlation function 



(46), which is related to the density-density correlation 



function ( 47 1 , displays a behavior characteristic of a liq- 
uid, we then conclude that the finite-momentum BEC of 
magnetic excitons is indeed a homogeneous phase. The 
disagreement between our results and the general analy- 
sis of Ref. ,35' might be related to the fact that here the 
electrons are restrict to the lowest Landau level. Recall 
that such a restriction, e.g., modifies the commutation 
relations between the electron and pseudospin density 



operators. 



VI. DISCUSSION 
A. Relation to previous work 

In this section, we briefly summarize some previous re- 
sults about the bilayer QHS and compare them with the 
ones derived here using the one and two-mode approxi- 
mations. 

ParlP31pj- p g ec i that the bilayer QHS at vt = 1 devel- 
ops a pseudospin spiral long range order at intermediate 
d/£ values. In this case, the main interlayer correlations 
are established between electrons and holes localized in 
different guiding centers [see Eq. (8), Ref. [H]. Inter- 
estingly, the excitation spectrum is gapped. It is argued 
that there is no fundamental reason for a Goldstone mode 
in this case (see note 14 in Ref. [2~Tj) . These two aspects 
above discussed suggest that the pseudospin spiral state 
bears some similarities with the finite-momentum BEC 
of bosons b. Moreover, it is also conjectured^ that the 
ground state is indeed given by a bound state between 
two pseudospin spirals with opposite winding direction. 
Recall that by splitting the finite-momentum condensate 
into two equal pieces, the total energy of the system de- 
creases, see Fig. |4j 

The ground state energy of finite size systems was cal- 
culated within the exact diagonalization technique. 20 It is 
shown that, regardless the size of the system, the ground 
state energy is almost constant for large d, a signature 
of the decoupling between the two layers. As we can see 
in Fig. HI the ground state energy of a finite-momentum 
BEC slowly varies for larger d/£. Moreover, the varia- 
tion decreases when the condensate is separated into two 
equal pieces. 

Nomura and Yoshioka also consider finite size systems 



and calculate the pair correlation function ( 46 ) via exact 
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diagonalization. 19 It is found that for d = 0.3 £ : both 
<7tt(r) and g^i(r) vanish as r —> but, for a larger 
d = 0.9 £, while g^^(0) vanishes, g-fi{0) is now finite. 
Concerning the large r region, both <7f-|-(r) and g^i(v) 
seems to be constant for d = 0.3 £ but, for d = 0.9^, 
they show small oscillations. Note that the pair corre- 
lation function, Fig. |11[ qualitatively displays the same 
features. 

Based on a Chern-Simons gauge theory, Bonesteel et 
alW* show that by approaching two composite fermion 
Fermi seas, there is always an instability towards the for- 
mation of composite fermion Cooper pairs. The theory is 
valid only in the large d region. The possibility of inter- 
layer composite fermion pairing is considered in Ref. [18] 
where some trial wave functions arc discussed. Assum- 
ing a p x — ip y pairing instability, it is shown that the two 
possible wave functions correspond to the (3,3,-1) and 
the so-called "strong" pairing (SP) states. The former 
phase displays a gapped (neutral) excitation spectrum. 
A qualitative phase diagram is also proposed and one of 
the possibilities is that the ground state changes as d/£ 
increases according to the following sequence: 111 - SP 
- (3,3,-1) state. Unfortunately, it is not clear how to 
compare a finite-momentum BEC of magnetic excitons 
with the SP and (3, 3, —1) states. 

Further support for pairing between interlayer compos- 
ite fermions is provided in Ref. 1161 However, the numer- 
ical results indicate that & p x + ip y pairing may occur 
instead of the p x — ip y considered by Kim et a/PS Mixed 
Fermi-Bose trial wave functions were then proposecP to 
describe the intermediate d ~ £ region, where the bosonic 
part is given by the 111 state while the fermionic one, by 
a paired composite fermion state. Such an approach in- 
deed follows the lines of an earlier work by Simon and 
coworkers pa where mixed wave functions are considered, 
but here the fermionic piece is given by two composite 
fermion Fermi seas. It is shown thalP the mixed wave 
functions with paired states provided a better descrip- 
tion for the intermediate d ~ £ region than the ones which 
do not include pairing. Again, it is difficult to compare 
this mixed Fermi-Bose wave functions with the finite- 
momentum BEC discussed here. We would like to point 
out that the description in terms of a finite-momentum 
BEC of magnetic excitons involves only bosonic degrees 
of freedom. 

Finally, studying the bilayer QHS within a Ginzburg- 
Landau theory, Ye and JiangS suggested that the ground 
state is given by a pseudospin density wave for d c \ < 
d < d C 2- In this case, the system undergoes two first- 
order quantum phase transitions: from the 111 state to 
a pseudospin density wave at d c i, and from the latter to 
two weakly coupled composite fermion Fermi liquids at 
d C 2 ■ The pseudospin density wave phase proposal is based 
on an earlier random phase approximation calculation,^ 
which finds that the (neutral) excitation spectrum has a 
minimum (magneto- roton) at finite momentum £q ~ 1.0 
and that the energy of this mode vanishes at d = 1.21. 
This feature indicates that a phase that spontaneously 



breaks translational symmetry may be realized. Such an 
inhomogeneous phase is studied, for instance, in Ref. [17] 
within the Hartree-Fock approximation. Interestingly, 
it is shown that the ground state energy is almost d/£ 
independent, similar to Fig. HI However, we should recall 
that the finite-momentum BEC is indeed a homogeneous 
phase, see Sec.[y] 



B. Consequences for the bilayer QHS at i>t = 1 



As discussed in Sec. Ill the two approximation 
schemes used throughout this paper, the one and two- 
mode approximations, impose some strong restrictions on 
the description of the finite-momentum BEC of magnetic 
excitons, but they allow us to carry out detailed calcu- 
lations in order to verify the stability of such a phase. 
Therefore, they should be seen as an initial approach 
to study the finite-momentum BEC. A more elaborated 
approximation scheme is needed. Since the two-mode 
approximation suggests that the system reduces its en- 
ergy by splitting the condensate into two equal pieces, 
a better approximation for the ground state should be 
a finite-momentum condensate such that all modes ±Q.; 
with Qi = Q and i = l,...n are (equally) macroscop- 
ically occupied, see, e.g., Refs. [55J for the case n —> oo. 
In particular, if n —¥ oo, the cylindrical symmetry of the 
quasiparticle dispersion relation f2 q would be restored. 
The implementation of such a scheme is rather involved 
and it will be deferred to a future publication. 

However, the results that we have derived so far al- 
low us to make the following statements about the bi- 
layer QHS: (a) There are strong indications that a finite 
momentum BEC phase is the most stable in the inter- 
mediate d/l region. Such a state bears a strong sim- 
ilarit y w ith the pseudospin spiral phase proposed by 
ParkPH (b) The instability of the zero-momentum BEC 
at d ~ 0.5 1, which we arrive at in Ref. HU indeed cor- 
responds to a first-order quantum phase transition from 
a zero-momentum BEC to a finite-momentum one. In 
principle, such a transition could be experimentally ob- 
served. 

It is also worth mentioning that (c) according to the 
two-mode approximation, a finite-momentum BEC with 
£Q = 1 is the lowest energy configuration only for 
d < 1.6 £, a value quite close to the critical d c where the 
incompressible-compressible quantum phase transition is 
experimentally observed. Curiously, the Fermi momen- 
tum kp of a composite fermion Fermi-liquid at v = 1/2 is 
£k F = 1. ESI The fact that £Q > £k F for d > 1.6 £ could be 
an indication of the incompressible-compressible phase 
transition. However, we should emphasize that this is 
just an interesting observation since, at the moment, it 
is not possible to identify a composite fermion Fermi- 
liquid phase within our bosonization formalism. 
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C. Open questions and next steps 

There are still a couple of open questions that we hope 
will be answered once the approximation discussed in the 
first paragraph of the previous section is carried out. ft 
remains to be verified whether: (i) The features found 
within the one and two-mode approximations [points (a) 
and (b) above] are robust, (ii) The momentum £Q as- 
sociated with the lowest energy configuration for a given 
d/£, Fig. pi either varies continuously with d/£ or is quan- 
tized, £Qo — 0,1,2,... (iii) It is possible to identify a 
finite-momentum BEC with £Q > 1 with two compos- 
ite fermion Fermi-liquids and/or with the pairing states 
proposed by Moller and coworkers.^ This would be an 
important step towards the determination of the nature 
of the incompressible-compressible phase transition. 

In addition to study a finite-momentum BEC where 
all modes Q^ with £Qi = £Q are macroscopically occu- 
pied, we also intend to consider the effects of a (small) 
finite electron interlayer tunneling, disorder (hopefully), 
and the electronic spin. Concerning the latter, there are 
some experimental evidence d 39 " 41 ! that the electron spin 
degree of freedom might be relevant for a complete de- 
scription of the bilayer QHS. For instance, it was recently 
reported 4 ^! that the critical d Cl where the incompressible- 
compressible phase transition takes place, increases and 
eventually saturates due to an increasing in-plane mag- 
netic field B||. In principle, the effects related to the 
electronic spin could be included in our analysis with the 
help of the generalized bosonization formalism^] which 
has been recently proposed by two of us to study the 
quantum Hall effect in graphene at v — and ±1. 

So far, we have focused on the elementary neutral exci- 
tations of the bilayer QHS. It remains to be checked how 
charged excitations are described within our bosoniza- 
tion approach. Such excitations are important when dis- 
order effects are taken into account (see, e.g., Refs. l43l 
and HU and the references therein). Two distinct cases 
should be considered: (a) Zero-momentum BEC of mag- 
netic excitons. Such a phase can also be seen as an XY 
pscudospin ferromagnet. — ' In this language, charged ex- 
citations correspond to topological (vortex) excitations 
called merons. There are four types of merons: with vor- 
ticity ±1 and electric charge ±e/2. For the single-layer 
QHS at v = 1, it was showrf 13 that a topological exci- 
tation (skyrmion) can be described as a boson coherent 
state. Due to the similarities between the single-layer 
and the bilayer QHSs, we expect that merons could be 
seen as a coherent state of bosons as well, (b) Finite- 
momentum BEC of magnetic excitons. Here the mapping 
into an XY pseudospin model no longer holds (see note 
14 in Ref. I2T]) and therefore, it is not yet clear whether 
charged excitations could also be described as a boson 
coherent state. 

Finally, concerning experiments, it would be interest- 
ing, e.g., to calculate the interlayer tunneling current 
for the finite-momentum BEC phase. The first theo- 
retical workg^ (clean limit) indicated that the bilayer 



QHS should display a Josephson-like effect, i.e., a zero- 
bias infinite tunneling conductance should be observed. 
Such a feature is related to the gapless linearly dispers- 
ing (neutral) excitation spectrum at low energies associ- 
ated with the Halperin 111 phase (zero-momentum BEC 
of magnetic excitons). Instead, an enhanced finite tun- 
neling conduc tanc e at zero-bias voltage was experimen- 
tally observed! 4 ^ In order to account for the experimen- 
tal features, disorder effects were then considered. At 
the moment, the experimental data have been under- 
stood within an XY pseudospin model with the tunnel- 
ing term (f2j) perturbatively treated and with disorder- 
induced merons phcnomenologically included in the elec- 
tron tunneling operator (for a review, see, e.g., Ref. l44l 
and the references therein). Interestingly, such a scheme 
indicates that by adding a parallel magnetic field Bn to 
the sample, the tunneling conductance versus bias volt- 
age data could provide a measurement of the gapless lin- 
early dispersing excitation spectrum. Again, the exper- 
imental data of Spielman et aZ.p who found some evi- 
dences for the existence of such collective excitation, were 
analyzed within the above XY pseudospin framework. 

In principle, our results are in contradiction with the 
experimental data^ since we have found a gapped phase 
in the intermediate dj £ region, where tunneling experi- 
ments were performed. However, note that according to 
our results, a true Josephson-like effect should occur only 
at very small d/£, where the zero-momentum BEC phase 
sets in. This is somehow in agreement with the exper- 
iments. Our next task is to verify whether the gapped 
excitation spectrum, Figs. |6| and [To} could account for the 
observed finite tunneling conductance at zero bias volt- 
age. In principle, we can calculate the interlayer tunnel- 
ing current (clean limit) within our bosonization formal- 
ism, treating the tunneling term (pi) nonpertubatively. 
Disorder effects could be included latter, for instance, 
following the lines of Ref. H3] In this way, we hope we 
can provide an alternative interpretation for the experi- 
mental data reported in Refs. |4] and HI 



VII. SUMMARY 

In this paper, we studied the bilayer QHS at Vr = 1 
within the bosonization methodp^l a formalism which al- 
lows us to properly treat the magnetic exciton as a boson, 
and we showed that the ground state of the system in the 
region d ~ £ (zero interlayer tunneling case) can be seen 
as a finite-momentum BEC of magnetic excitons. Our 
findings are in agreement with previous results which sug- 
gest that an intermediate phase may show up between 
the (incompressible) Halperin 111 state (ground state 
for small d/£) and the (compressible) composite fermion 
Fermi-liquids (ground state for larger d/£). 

The stability of the finite-momentum BEC has been 
analyzed via two distinct approximation schemes: the 
one-mode approximation, where it is assumed that the 
bosons macroscopically occupied only one momentum 
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Q = Qx with £Q ^ 0, and the two-mode approximation, 
where both ±Q modes with Q = Qx are macroscopically 
occupied. We have found that such a phase can be real- 
ized as long as the excitation spectrum is gapped. The 
comparison between the ground state energy curves in 
terms of the ratio d/£ for configurations with different 
£Q as well as the analysis of the quasiparticle excitation 
spectra provide strong evidences for a first-order quan- 
tum phase transition at small d/£, i.e., a transition from 
a zero-momentum BEC, a phase that we had already an- 
alyzed in Ref . IT41 and that corresponds to Halperin 111 
state, to a finite-momentum BEC. We hope that such 
first-order quantum phase transition can be experimen- 
tally observed in the near future. 

As a final remark, we would like to emphasize that 
the bosonization method introduced in Ref. T3] can be 
used to study both the single and double-layer QHSs at 
v = 1. In other words, we can describe both systems 
using the same degree of freedom, the magnetic exciton, 
in the limit where this object can be treated as a boson. 
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Appendix A: About the bosonization scheme 

In Ref. 1131 it is shown that a bosonization formalism 
for the two-dimensional electron gas under a strong mag- 
netic field (single-layer QHS at v = 1) can be developed 
following the lines of the bosonization scheme used to de- 
scribe one-dimensional electronic systems. It is interest- 
ing that this formalism also gives quite reasonable results 
for the bilayer QHS at vt = 1 even though such a system 
is in a limit very far from the one considered in Ref. 1131 
In this section, we start providing some heuristic argu- 
ments which tell us why the bosonization scheme^ can 
also be employed to study the bilayer QHS at vt = 1 
and later, we show a simple calculation which corrobo- 
rates such arguments. 

The bosonization method for the single-layer QHS at 
v = 1 is based on the following points: The ground state 
of the (noninteracting) two-dimensional electron gas at 
v = 1 is the quantum Hall ferromagnet |FM), Fig. [2] (a), 
the reference state. Elementary (neutral) excitations are 
electron-hole pairs or spin flips, Fig. [21(b), which can be 
obtained by applying the spin density operator S~ into 
|FM), i.e., |*) ex 5"|FM). Although the commutation 
relation between the spin density operators S + and S~ 



(projected into the lowest Landau level), 



Ivq ' ^k 



J 2 qk*/2 



Pt 



(q + k)- e ' 2 ^/ 2 pi(q + k) (Al) 



with q = q x + iq v , differs from the canonical commu- 
tation relation between a creation and an annihilation 
boson operators, it is still possible to define boson oper- 
ators b as in Eq. |9| if we follow the lines of Tomonaga's 
procedure's for one-dimensional electron systems. Using 
the fact that 



p r (q)«(FM|p r (q)|FM)=iV$ q ,o, 



^(q)«(FM|p 4 (q)|FM)=0, 



we realize that Eq. ( Al ) assumes the form 



Pq.Sk] ~ A^expH^/^q+k.O, 



(A2) 



(A3) 



which now resembles a canonical commutation relation 
for boson operators. In other words, as long as we are 
close to the |FM) state, i.e., the number of bosons in 
the system is small, electron-hole excitations (magnetic 
excitons) can be approximately treated as bosons. 

Turning to the bilayer QHS at vt = 1, we note that 
such a system is very far from the limit discussed above 
because the configuration vt = v\ + v\, = 1/2 + 1/2 = 1 
corresponds to a system with N^/2 bosons, see Fig. fl] 
(b). As we will see below, this is indeed a very special 
configuration where density fluctuations guarantee that 
the relation (A3) still holds. 

Let us consider as a reference state the zero- 
momentum BEC of magnetic excitons illustrated in Fig.fTI 
(b) and identify the electronic spin degree of freedom of 
the single-layer with the pseudospin one of the bilayer 
QHS. In this case, we have 



Pt(q) ~ (pt(q)) +fyt(q), 
Mq)~ W(q)) + %(q)- 



(A4) 



Since (p^(q)) = (pj,(q)) = ( N^/2)6 Qt o, we note that here 
the equivalent of Eq. (A3) vanishes, and therefore the 
bosons b are no longer well defined. However, since the 
total filling factor is fixed, we have 



<W(q) = -%(q) = fy(q)> 



(A5) 



which indicates that a finite (approximate) commutation 
relation can still be obtained if we now consider density 
fluctuations in Eq. ( A4 ) . Due to the relation between the 



magnetic exciton momentum q and the guiding centers 
associate with the electron and the hole, we have 



Sp(q) = <5/9local(q) + <5pnonlocal(q) 



(A6) 



where <5pi oca i(q) and 6p n oniocai(q) correspond to density 
fluctuations on the same and different guiding centers 
respectively. Moreover, we can also write 



^Piocai(q) ~ PoSq,o, 



(A7) 



15 



with po being a constant, which implies that the relation 
(A3) still holds for the bilayer QHS at vt = 1- 



In order to see that (A6) and (A7) are indeed quite 



reasonable assumptions, let us calculate the density fluc- 
tuations 



S Pa (q) = V(P 2 M)) (Pa(q)> 2 



(A8) 



within the two-mode approximation. It can be seen as a 
kind of self-consistent check of the above assumptions. 

In the Bogoliubov approximation, the Fourier trans- 
form of the a-electron density operator reads [see 
Eqs. (25) and (26) from Ref. [T3] 

p Q (q) - <W<*,,oJV» " ae~ ( ^ )2/4 [e- MqAQ/2 0Vb&^ +q 
+ e iQqAQ / 2 A /^o6 t _ Q+q 



\p -mqA(Q+k)/2,t u 

/ / e °Q+k+q°Q+k 

k^0,-2Q 



(A9) 



Since 



1/2 



ri/2 



(bl Q+cl )=5 clfi N 1 a /2 + S (ii2Q N 1 Q /2 1 



and 



( b Q+k+q 6 Q+k) = <5q,0(^ +k 6 Q +k) +5 q ,2Q(63 Q+k 6Q+k) 
+ <J q ,-2Q<6L Q+k 6QH-k), 

where the expectation value is taken with respect to the 
ground state of the bilayer QHS, i.e., the vacuum for the 
bosons a, see Eq. (36), it follows that 



(Pa(q)) - <5 q ,o^ - (<5 q , 2Q + 5 q ,_ 2Q )e-^ 2 F Q (q). 



(A10) 



Here, 



_l_ y^ r „,q i- - 

k^0,2aQ 



('■'-■ -|Vi(k)i) 2 (k) + u 3 (k)w4(k)] 



with Uj(q) being the Bogoliubov coefficients, Eq. (C5). 
Concerning (p 2 (q)), after some algebra, it is possible to 
show that 



(pa(q)Pa(q)> = <5q.O 



1 



r^J + Po 



3q,2Q + °q,-2Qj 



x e ^Q) 2 [F Q 2 (q)+ Pl (q)], (All) 



where 

(°o = J2 K ( p ) + W 3(P)] 2 + [ui(p)«a(p) 



+u 3 (p)u 4 (p)] 2 + 2 Mp)v 2 (p) + «3(p)w4(p)] 2 



+2 [«i(p)vi(p) + W3(p)w 3 (p)] 2 , 

pi(q) = Yl cos ( q A p ) MpWp) + w 3 (p)w 4 (p)] 2 

p 

+ Mp)v 2 (p) + u 3 (p)vi{p)} 2 . 

Therefore, Eqs. {A5}, |A8| , ( |ATo| , and ( |A"TT) ) yield 



^P 2 (q) = <5q,0P0 + [^q,2Q + <5 q -2Q] 



-2(£Q) 2 



Pi(q)- (Ai2) 



Note that the first term of the above equation can be 



identified with Spi oca i(q) in (A6) while the second one, 



with i5p non i oca i(q). Moreover, one can easily see that 
pa > e _2 ( £< 3) pi(±Q), which corroborates the fact that 
<5pnoniocai(q) can be neglected with resp ect to <5pi oca ^(q). 
In the one-mode approximation, Eqs. ( A10) and (All ) 
reduce to 

1 



(Pa(q)) =Sq, -N^, 



(pa(q)p<r(q)) = 5 q , 



1 !\r2 , ! V £ k( £ k - ^k) 

4^> + 2 ^ 

k#Q 



n£ 



Although ^Pnoniocai(q) = 0, the bosonization scheme^ 
still holds because the relevant term 5pi oca i(q) is finite 
within this approximation. 



Appendix B: Alternative effective boson models for 
the bilayer QHS 

In this section, we briefly comment on some different 
effective boson models proposed to describe the bilayer 
QHS. 

In Ref. Hll the effects of the electron spin degree 
of freedom are taken into account. Here the original 
fermion model is mapped into an effective lattice spin- 
pseudospin model, which is then analyzed within a gen- 
eralized Schwinger boson mean-field theory. Finite tem- 
perature properties are calculated, for instance, the tem- 
perature dependence of the spin and in-plane pseudospin 
magnetizations. A proper comparison between our re- 
sults and the ones of Ref. |47] will be possible only af- 
ter including the electronic spin in our formalism, see 
Sec. |VIC| We would like to recall that our approach 
is based on a direct mapping of the interacting fermion 
model (fTl) into the boson model (12), no lattice degrees 
of freedom are introduced. 

Tieleman and co-workers-^ derived an effective boson 
model for the bilayer QHS also following the ideas of the 
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bosonization schemep^ but with some important differ- 
ences: The single particle electron states considered are 
no longer the pseudospin up and down lowest Landau lev- 
els, see Fig. Ill but symmetric and antisymmetric linear 
combination of these states. Instead of defining the boson 
operators with respect to the quantum Hall fcrromagnet 
|FM), see Sec. II A the reference state is the completely 



filled symmetric state |SYM). Therefore, the boson op- 
erators introduced in Ref. EH hereafter called (3, differ 
from the bosons b discussed in Sec. |II A| Most impor- 
tantly, the ground state corresponds to an almost filled 
symmetric state instead of the BEC of magnetic excitons 
considered within our formalism. Since the procedure 
adopted in Ref. [48] to derive an effective boson model 
is not fully consistent with the bosonization scheme,^ 
below we briefly revisit its derivation. 

Formally, the bosonization schemes of Refs. fT3l and l48l 
are identical. Therefore, the bosonic representations of 
the electron density and pseudospin density operators are 
given by Eqs. (27)-(29), and (31) of Ref. 03] with the 
replacement b — > (3. Substituting these expressions in 
Eq. (7) of Ref. SSI we arrive at 



jjSAS 



H 2 + H 4 + Hq , 



(Bl) 



where #2,4,6 respectively correspond to terms with 2, 4, 
and 6 boson operators /3. In particular, 

#2 = ]T e q /3 q /3 q + I (A q /3t/3l q + A q /3_ q /3 q ) , (B2) 



with 



= t + A q -^< 



-{lkf/2 



v c (k) 



+ J2 2fo(fc)e~ ( ' fc)2/2 sin 2 (k A q/2), 

k 

A q = N^e-^ 2 / 2 v c (q); (B3) 



A q = A q - Y^ Vc(k)e-W I 2 cos(k A q). 



Here t and v / c (k) respectively correspond to Asas and 
v o/z(k), see Sec. [TTl Comparing Eqs. (B1|-(B3| with 
equations (19)-(21j of Ref. |48j one can see that the 
former have extra terms. This is related to the fact 
that Tieleman et al. approximated the bosonic repre- 
sentation of the operator ^(k) by linear terms, i.e., 
S x (k) ~ pi + /3_kj while here the complete bosonic rep- 
resentation, which in addition includes a cubic term in 
boson operators /?, is considered. We should mention 
that the presence of such a cubic term in the bosonic rep- 
resentation of S x (k) is important because it guarantees 
that the bosonic representation of the electron density 
and spin density operators obey the correct commuta- 
tion relations, the so-called lowest Landau level algebra, 
see Sec. II. D of Ref. 13 for details. In order to perform a 
proper map of the original fermion model into the boson 



one, we should consider the complete bosonic expression 
of the electron density and spin density operators. After 
that, approximations can be employed. 

Therefore, if we define boson operators with respect 
to the state |SYM) and try to be consistent with 
the bosonization formalism,^ we then obtain a non- 
Hermitian effective boson model to describe the bilayer 
QHS. Recall that the procedure adopted in this paper 



yields a Hermitian boson model, Eq. (12 1 



Although the Hamiltonian (Bl) is non-Hermitian, let 



us for the moment assume that this is only an artefact 
of the bosonization scheme and determine the spectrum 
of the elementary excitations J7 q within the harmonic 
approx ima tion, i.e., H§ AS m #2- It is possible to diago- 
nalize (B2) as done, for instance, in Ref. 49. One finds 



ft q = 



v, A a A, 



q' v q- 



However, f2 q is not well-defined for small momenta: it is 
easy to see that e q — A q A q < for q — > 0. One concludes 



that the restriction of Hg AS to the quadratic term ( B2 ) 



is not a good approximation. A proper analysis of ffg^ 
should take into account the H4 and Hq terms. Finally, 
one should mention that the approach of Ref. 08j is more 
suitable for larger AgAs values while our formalism, for 
the opposite limit. 

Interestingly, apart fro m t he coefficient of the /3 q /3l q 



term, the Hamiltonian (Bl ) corresponds to the effec- 
tive boson model proposed in Ref. 50: Here the en- 
ergy of the elementary excitations is given by J7 q = 
(e„ — A 2 ) 1 / 2 , which agrees with the diagrammatic cal- 
culations of Fertig:-3I Recall that, in this case, the spec- 
trum has a magnetoroton minimum which vanishes for 
ds 1.2^ (t = case). 



Appendix C: Details: two-mode approximation 

In this appendix, we provide the full expressions of 



some quantities which appear in Sec. 



After the substitution (|27h, Eq. (15) acquires the form 



IV 



K = K ° + \ E' [<( fe Q +q 6 Q+ q + & -Q-q &t -Q- q ) 



£ q ( fo Q-q fo Q-q + & -Q+ q & -Q+q) 



7qOQ +q k-Q+q + h -C- + k-Q-qfrq-q + k.C.) 



A q (Oq_ q 



h.c.) 



A q( &t -Q + q^Q-q+ h - C -) 



t »,t 



£ q (& T _ Q+q ^_ q + h.c. + 6 T Q+q 6 T _ Q _ q + h.c 



(CI) 
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The restriction on the sum over momenta indicates that 
the modes that satisfy the condition ±Q ± q 7^ ±Q are 
not included. Here 



with i = l,2, 3, 4, ft q = Q+ (i = 1,3), fi" (i = 2,4), 
and on = 1 (i= 1,2), -1 (i = 3,4). 

Appendix D: Pair correlation function 



K = 2N NqV2q{Q, Q) + (wq — ^)(N + N ) V^ e q In this section, we present the full expressions for the 



7q = WN N [« q (-Q, Q) + « 2Q (q, Q)] , 
A q = 4A^ q (Q,Q), A q = 4/VoU q (Q,Q), 

e q = 4v/iVoiVo Uq (-Q,Q), (C2) 



and e q are defined in Eq. (18). Setting N — N and 



using the fact that tiq = Nq/Nb = 4tt1 2 Nq, we can write 

P 2 e -(lqf/2 



A, 



-(2no) 



{<; 



-[l-e-^ d ^cos(qAQ)], 



2 -(lqf/2 

£ q = ^-(2n )^— [cos(q A Q) - e-<W], (C3) 



e 2 e -(iqf/2 



-2{lq)d/h 



cos(q A Q), 



which are useful expressions for the numerical calcula- 
tions. 



Concerning the chemical potential, Eq. (41 ), the terms 
/j,q and /xi are defined as 



Mo 



— E 



1 1 

1 i a 1 Let 



[A q (e q -A q ) + 7 ^-^] 



- 2 ^ + ^(^-^.)[^'V^-A q , 



+ (7q£q ~ 7qCq)(7 q e q + 7qA q - 2A q £ q )] , (C4) 



Ml = 2-/V" U 2 q(Q,Q) = ^o«0 



e- 1 e" 2 W) 



d 2 ZQ 



1-e- 



-2(lQ)a 



Finally, following the procedure described in Ref. [3T| 
(see Chap. 3) and after some algebra, it is possible to 
show that the Bogoliubov coefficients Ui(q) and Uj(q), 



the elements of the 4x4 matrix M q , Eq. (37), are given 
by 



«?(q),o?(q) 



1 
4^ 



(e, + a»A q ± Q q ) 



(-1) 1 

" — [7q (7q£q ~ Aq£q) - a 4 A q A q 



4^q^q 



A q e q (a,e q + A q ±a 4 r! q )] (C5) 



pair correlation function ( 46 1 , calculated both at the one 



and two-mode approximations. 

In order to calculate the static structure factor ( |47| , 
it is necessary to project the product p{— q)p(q) into the 
lowest Landau level. We have [see Eqs. (4.17) and (4.18) 
from Ref. [36] 



p(-q)p(q) = K-q)p(q) + JV[1 - exp(-|Zq| 2 /2)], (Dl) 

and therefore Eq. ( [47] ) assumes the form 

5(q) = 5(q) + [1 - cxp(-|/q| 2 /2)] (D2) 

with 5(q) = (l/iV)(p(-q)p(q)) - iV£ q , being the pro- 
jected static structure factor. From Eq. (10), it follows 
that (q ^ 0) 



S(q) 



4e -|<9| 2 /2 

2^ sin(q A p/2) sin(q A k/2) 



N, h 



P .k 



:<4 



k-q & k^ +q &p)- 



(D3) 



Here the exp(— |/g| 2 /4) factor is restored in the bosonic 
representation of the electron density operator, see 
Eq. (27) from Ref. M 

Using Wick's theorem, the expectation value in the 
above equation can be easily calculated. P erfo rming the 
replacement 6q , 6q 
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and keep- 
it is pos- 



'No as done in Sec. 
ing only the terms with two boson operators 6', 
sible to show that in the one-mode approximation the 
projected static structure factor assumes the form 

5(q) « 27i e-l / «l 2 / 2 sin 2 (qAQ/2)(t; 2 + U 2 -2 Uq t, q ) (D4) 

with the Bogoliubov coefficients u q and w q given by 
Eq. (23). Note that 5(q) vanishes if Q = implying 



that, within this level of approximation, 
g(r) = 1 - cxp[-r 2 /(2Z 2 )] 



(D5) 



regardless the value of d/L Eq. (D5) is nothing else but 



the pair correlation function for the single-layer QHS at 
v = \. 

The same procedure can be used to calculate g(r) in 
the two-mode approximation. After some algebra, we 
find that 



5(q) 



2noe-l^ 2 / 2 sin 2 (qAQ)[ 

[?ii(q) + «i(q) - u 2 (q) - v 2 (q L )} 2 

[«a(q) + « 3 (q) " «4(q) - «4(q)] 2 ] , (D6) 



where the Bogoliubov coefficients u.;(q) and Vj(q) are now 
given by Eq. (IC5 1 . 
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